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Abstract. Asymmetric incommensurable torque signals can be achieved by exploiting magnetic 
nonpolar repulsion techniques, see also [1], aiming at designing kinetic energy harvester (KEH) 
systems. This paper confirms by using three different simulation methods that an asymmetric 
incommensurable torque signal is feasible to be reached, exploiting non-conservative magnet 
fields. All shown calculations are in fully agreement with the classical Electromagnetic Theory. 
The first method is based on time consuming full 3D FE generated torque and force signals. The 
second method is using also FE generated torque and force signals, but they are calculated only 
in orthogonal directions, approximated with Fourier Series and normalized to provide 2D 
stiffness signals. The third approach focuses on filamentary current rings to calculate the 
magnetic field of permanent magnets (PM) to generate in a purely analytically way the same 
force and torque signals as in the previous methods. However, only speculation about the energy 
source as required to maintain the forecasted persistent torque can be given. 

1. Introduction 

Energy Harvesting is a technology for capturing non-electrical energy from ambient energy sources, 
converting it into electrical energy and storing it to power wireless electronic devices [2], [3], [4], [5]. 
The process of capturing mechanical energy such as shocks and vibrations is a subject area of energy 
harvesting requiring specific types of devices, so called kinetic energy harvesters (KEH) and there are 
many types of KEH’s [5]. For the devices described in this paper we employ still the naming convention 
KEH, as effectively rotor mass and rotor inertia are used to start an everlasting oscillation in theory. 
Microscopically, there might be yet unknown sources why such a persistent motion is possible to obtain 
[6]. In chapter 2 asymmetric incommensurable torque signals are modeled using the three different 
approaches. In chapter 3 a proposition for a KEH prototype buildup is discussed. In the conclusions, 
chapter 4 we summarize the findings. 

2. Design of asymmetric incommensurable torque signals 

The force or torque of two permanent magnets (PM) can be calculated using Lorentz force law 

Fmag = </(E + (vx B)) = q(vx B) (1) 


where F mag is the resulting force, E the electric field-, B the magnetic field- and v the point charge q 
velocity trajectory-vector. As in these KEH systems, velocities are extremely low in respect to the speed 
of light, E can be disregarded (which is always assumed in calculations for electrical machines, see for 
instance [7]). Generally, magnetic forces will do no work on isolated electric charges [8], since 


1 / 1 / „ = 
vv mag 


mag 


dl 


d\=vodt 
- > 


j q(y x B) ■ v 0 dt = 0 



The result in (2) follows when v and the created path velocity vector v 0 are equal, which is always the 
case using a point charge. However, for the PM-based KEH systems proposed in this paper, we show 
that the calculated total magnetic work W mag is unequal zero by applying three different methods and 
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using (a) well-deliberated closed trajectories T, and (b) instead of a point charge, a disk PM or a current 
loop (both having geometrical extensions). So, in contrast to (2), 

W m ag = Fmag m dl =1= 0 (3) 

The simplest possible configuration of a KEH system that can exert asymmetric incommensurable 
torque is shown in Figure la. Figure lb, c show the normalized 2D stiffness for the torque- and force- 
signal assuming an airgap of 1mm between stator-rotor PM set. The generation of such amplitude 
diagrams can for instance be approximated with Fourier Series. It is sufficient to consider only terms up 
to the forth order, see also [1]. These normalized signals for torque and force components, given by 

Trad (0, Z) = C r fz (0, (4) 

4 x (0,z) = 4/ F (0,z) (5) 

respectively, are stiffness functions, where the torque component C r ( Nm/rad ) and its counterpart force 
component k r (N/m) can be easily measured as the maximal break-torque respectively break-force. 
Introducing a harmonic axial cam movement of the rotor as 

z F (0) = A co sin(a) co 0 + <f co ) (6) 

where the lateral movement z R is constrained to the rotating DoF 0 with a given cam amplitude A co , an 
initial offset <f co and an angular velocity a) co . 



20 Magnetic Spring Function • Fore* Component 


20 Magnetic Spring Function - Torque Component 


Rotor disk 
magnet with 
<Z>d R ,h R (2) 
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Stator disk 
magnet with 

Wji h s ( 1 ) 




(a) (b) (c) 

Figure 1. Concept view (a) of rotary-translatory nonlinear PM spring system, with a rotary (0), an axial (z) and 
a radial (r) DoF. Normalized 2D torque spring and force spring amplitude for 0 and z direction shown in (b) 
resp. (c). Shown data represents a disk PM with diameter 10mm, height 5mm, residual mag. N52 (B = 1.4457). 




In the following simulations, the values in (6) are kept as a) co = 3 and <f co = 0 , to make the case 
simple. The setup is illustrated in Figure 2, showing also a cam function with a realistically small cam 
amplitude of A co = 1 mm and an asymmetric rotor to stator displacement of 9 mm. 

Developed View of PM Stator and Rotor Distribution (axR3S3sl) so 
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Figure 2. Concept buildup of a rotating and axial wobbling 3PM rotor set and an asymmetrically placed 3PM 
stator set in a single line distribution (axR3S3sl) with developed view (a) and front view (b). 

Figure 3a shows the optimization process for creating an unbalanced torque signal by sweeping 
the stator-rotor asymmetry from d asym — 0 ... 20 mm. Keeping the stator and rotor PM in the same 

z-axis with an enabled cam movement at 0mm, no asymmetric torque nor force signal in one 
revolution occurs. As soon as an asymmetry (z offset) between rotor and stator PM is introduced, the 
(normalized) torque signal will increase, and at half the PM-diameter it reaches a maximum, followed 
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by a decline. When offset reaches 15mm, no torque signal at all is generated, due to the lacking mutual 
PM magnet influence. The torque signal (blue) is depicted in %, showing 50% as a balanced torque 
signal. No energy surplus is occurring, e.g. the created force signal needs accordingly to be damped 
to generate an energy surplus (see also chapter 3). Figure 3c shows generated torque signals with FE 
simulation (red) and Fourier Series approximation (blue) plus analytical model using Lorentz 
approach [9], [10], [11]. In all three cases, a net torque of all summed up rotor disk PMs is clearly 
positive (61%) over one full revolution; peak torque is reached at 10°, 130° and 250°. 



Figure 3. Diagram (a) depicts a stator PM distance sweep from 0...20mm; resulting torque signals in % (left 
axis) and normalized (right axis). FE model in (b) shows the PM setup for the full 3D FE simulation; (c) shows 
the corresponding created torque signal of all summed up 3 rotor disk PMs applying the 3 different methods. 


This behavior is most intriguing, as by using solely well accepted classic EM theory such as Lorentz 
force [9], Biot-Savart or Mutual Inductances [10] approach, an energy surplus is clearly forecasted. An 
energy source must therefore be present somehow, and some conjecture about such a source is given in 
[1] and [11]. 

3. Proposed system for buildup 

A system concept buildup is proposed in Figure 4, with focus of keeping the cam amplitude low (with 
A co < 3 mm) and allowing an angular cam velocity coupling of co co = 1 only. 


Developed View of PM Stator and Rotor Distribution (axR2S3ml) 
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Figure 4. Concept buildup with developed view (a) and front view (b) of the proposed arrangement with two 
rotor PMs (green and numbered) and three stator PMs (blue) using multi line PM distribution (axR2S3ml). 


Using the PM configuration and the sketched cam movement in Figure 4, the motion profile (torque 
DE) with constraining force created from lateral movement (6) follows (using Lagrange approach): 

11 1 

70" + -A C0 2 ma) C0 2 <p" + —A C0 2 m a) co 2 cos(2 ai co 0) 0" - - A C0 2 ma ) C0 3 sin(2 a) C0 cp) (j)' 2 + r D (0', 0) + r E (0) = 0 (7) 


In case the torque signals from the constraining force are neglected from (7), e.g. A co 
DE can be written much simpler as 

/0" +T D (0',0) +t s (0) = 0 


0, the torque 



with J denoting the inertia, t d the torque due to friction and the asymmetric incommensurable torque. 

The simulation model confirms that, with a realistic rotor mass of m — 200 g (J = 45 figm 2 ), a 
persistent motion is generated, as shown in Figure 5d, having applied a nonlinear friction model 
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(compare Figure 5c, used also in [1]), and an approximated SDoF dynamic model for which the 
acceleration and deceleration due to the lateral cam movement is taken into account (solid lines in Figure 
5b). Also, the nonlinear torque signal r D should be considered as 

Tfl(0',0) = T Dradi 0') + *Dax(<P'' </>) = (i D r adW>') + D ax {(f)')\ f F (<p, Z R ($)) \)(f)' (9) 


with torque friction component D rad ((p f ) and Z) ax (0'). Both components are described with a nonlinear 
friction model, depicted in Figure 5c. In (6), D rad (<fi ') 0' is the torque friction for revolving the rotor 
(which could have been also modeled with a linear function, as the rotor is revolving at steady state well 
above 0' > 10 rad/s), and D ax (<p r ) | / F (0, z F (0))| 0' is the additional torque friction introduced 
when rotor PMs are in proximity of the stator PMs. The normalized force / F (0, z F (0)) has been 
introduced earlier (see Figure lc, as we deal with friction, the absolute value of f F must be taken, as the 
rotor exerts on either cam well a positive damping force resp. damping torque; resulting constraining 
force F ax depicted in Figure 5a). This constraining force cannot be eliminated (but it could be minimized 
with an elaborate stator-rotor PM arrangement), as generated forces in lateral z-direction for each stator- 
rotor PM rendez-vous are changing sign (compare also Figure 5a and Figure 4, PM rendez-vous at la 
and 2a are doubling negative z-forces and after 180°, at lb and 2b, these forces changing sign). 

For creating an everlasting motion it is important to note that the friction torque signal over one full 
revolution must be smaller than the asymmetric incommensurable torque signal created by the stator- 
rotor PM pair arrangement (shown in yellow, Figure 5a). The magnetic energy W mag of this dynamic 
system in steady-state for one revolution must fulfil 

r2n r2n 

W mag =\ ^(4>)d(p = T D (4>',4>)d<p' (10) 

Jo Jo 

The complexity of this proposed configuration is modest and expenditures for prototyping is small. 
Even though the lateral constraining z-forces are symmetrical, this lateral energy component is not 
canceled out over one revolution and a low friction cam system must be used to guarantee that the energy 
component from the asymmetric incommensurable torque signal is larger. 
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Figure 5. Diagram (a) depicts the resulting normalized torque and force signals from proposed setup in Figure 4, 
using A co = 3mm and a) co = 1 (<f co = 0) applying the 3 different calculation methods; (b) shows motion functions 
of approximated cam signal with and without ( A co = 0) constraining force (and keeping r z zero); (c) friction 
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model from [12] for radial and axial friction (in set rotor position without additional lateral load); (d) resulting 
motion signals. 

The torque asymmetry can be further increased using different PM shapes (for instance, using block 
magnets, an asymmetry of nearly 60% could be achieved and larger force signals will result). In addition, 
such a system can be easily stacked in z-direction, or the coupling increased to a) co = 2 (making it 
necessary to double the PM arrangement frequency on the surface for stator and rotor, compare Figure 
4) to basically unfold an energy surplus as desired. 

4. Conclusions 

It is stated that, in theory, by combining trajectories of a rotor with accordingly stator-rotor PM 
distributions, an asymmetric incommensurable torque signal is feasible. The presented trajectories have 
been limited only in rotational- (0) and translational- (z) direction, excluding movements that involve 
also radial direction (changing rotor diameter distance r, see also Figure la). However, we showed with 
three different approaches that the same asymmetric incommensurable torque signal can be numerically 
simulated based on well accepted classic EM theory only. While keeping a cam harmonic, the generated 
force signal must be damped in all cases - by using a low friction cam setup, whether using a PM setup 
from Figure 2 or from Figure 4 for creating a persistent motion. Also important in this study is the 
foresight, that magnet field distributions can exert energy and are not per se conservative fields. 
According to the presented results, the magnetic field itself is neither conservative nor non-conservative; 
it all depends on the closed path trajectories followed through the field. Non-conservative trajectories 
can be generated which allow netto work while rotating a rotor over one full revolution. The fundaments 
for the forecasted energy source, which is necessary to deliver this netto work, remain a research 
question, and further open-source experiments are indispensable to verify the presented theoretical 
claims. 
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